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1. INTRODUCTION 

Over the course of centuries there has been enormous 
effort devoted to theories and experiments for classical 
systems in external gravitational fields [J [1] 0 0 • Two 
recent experiments purport to give evidence for a grav¬ 
itomagnetic field, which is the gravitational analog of a 
magnetic field due to the motion of source matter: the 
LAGEOS experiments detect gravitomagnetic effects on 
the orbital motion of two earth satellites, and the Grav¬ 
ity Probe B experiment detects gravitomagnetic effects 
on the precession of four gyroscopes in earth orbit 0 0. 

Much attention has also been given to the theory 
of quantum fields in classical background gravitational 
fields, in particular regarding Hawking radiation by black 
holes 0. Notable work has been done on quantum sys¬ 
tems in classical background gravitational fields, for ex¬ 
ample on neutrons in the earth’s field 00, and more 
recently for atomic beam interferometry; this last work 
is largely based on semi-classical calculations of phase 
shifts along the classical atomic trajectories 0- Much 
less work has been done on gravitomagnetic effects on 
fully quantum mechanical systems. Experiments to de¬ 
tect such effects should clearly be expected to be quite 
difficult, but would be of fundamental interest (llj[ia|. 

In this work we give a systematic treatment of a scalar 
or spin zero quantum particle in a combined electromag¬ 
netic and weak gravitational field. We describe the parti¬ 
cle with the generally covariant Klein-Gordon equation, 
minimally coupled to the electromagnetic field in stan¬ 
dard fashion [10. The weak gravitational field is nat¬ 
urally treated according to linearized general relativity, 
and we also limit ourselves to slowly moving sources. In 
this case the gravitational field equations are quite analo¬ 
gous to those of classical electromagnetism [14j ; we refer 
to this as the gravitoelectromagnetic or GEM limit. 

We proceed by first casting the Klein-Gordon equation 


in an exact Schroedinger equation form (SEF), which cu¬ 
riously does not seem to appear in the literature jl5j . 
The SEF of course lends itself well to considerations of 
the non-relativistic limit for the quantum particle, which 
we then study to leading order in the relevant energies 
0| . The reader not interested in the theoretical devel¬ 
opment may skip directly to the results in section 5, in 
particular Eq. (5.1). 

Our main focus is on gravitomagnetic effects on a quan¬ 
tum system, so we briefly consider a number of possibly 
observable physical effects, in particular using a fictitious 
spin zero electron bound to an atom in a gravitomagnetic 
field. (That is we ignore spin.) It is conceivable that some 
novel interactions we derive, involving the product of the 
EM vector potential and the gravitomagnetic field, might 
be of interest. In considering such physical effects we only 
present parametrizations and rough numerical estimates 
since it is beyond the scope of the present work to an¬ 
alyze specific experiments; we leave such considerations 
as a challenge to experimentalists. 

In a subsequent paper we will consider spin 1/2 parti¬ 
cles and compare the results to those of the present work, 
especially those related to gravitomagnetic interactions 
with spin angular momentum. 

It should be emphasized that here we treat gravity en¬ 
tirely classically, so our work does not relate to quantum 
gravity or quantum spacetime 00 - 


2. GRAVITY IN THE 

GRAVITOELECTROMAGNETIC (GEM) LIMIT 

In this section we review very briefly linearized general 
relativity for slowly moving matter sources such as the 
earth, paying special attention to the gravitomagnetic 
part of the field 10 - Note that it is the gravitational 
source matter that we assume is moving slowly, not the 
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quantum objects in the field. We use h = c = 1, but 
retain dimensions for G. 

For a weak gravitational field the metric may be writ¬ 
ten in Lorentz coordinates as 

= r Uiv + fytuo (2-1) 

where is the Lorentz metric and h is a small per¬ 
turbation. For ordinary matter such as that of the earth 
the energy momentum tensor is well approximated by the 
matter or “dust” tensor 0, 

T^jy = pU^Ui!. ( 2 - 2 ) 

Here p is the scalar density of the source matter and u M 
is its 4-velocity. To first order in h^ the field equations 
are 0 


<9 2 (/i m „ - -Tjftuh) = -167 xGpu^Uv. (2.3) 

Here d 2 = d 2 — V 2 is the d’Alembertian operator and 
we have used coordinate freedom to impose the Lorentz 
condition 


{h,j,v ~ 2 = 0. 


(2.4) 


The single slash denotes an ordinary derivative. 

For slowly moving source matter we may ignore second 
and higher order terms in the source velocity. Then the 
linearized field Eq.s (2.3) imply a simple form for 
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h k = h 0k ,^ = h 00 /2, (2.5) 


where 0 is the Newtonian or gravitoelectric potential, and 
h k is the gravitomagnetic potential. The fields 0 and h k 
obey 


—<9 2 0 = AnGp, — d' 2 h 3 = —167 rGpv 3 , (2.6a) 

40- V-h = 0, h = 0. (2.6b) 


Equations (2.6) are of course almost identical to the 
equations of classical electromagnetism (EM) for the 
Coulomb potential <f> c and the vector potential A 3 , 

—<9 2 0 c = 4irp qi — d 2 A 3 = AirpqV 3 , <j> c + V ■ A = 0. 

(2.9) 

Here p q is the charge density and the Lorentz gauge is 
imposed. The only difference between EM and GEM 
equations is a factor of —4 associated with the gravito¬ 
magnetic potential. Moreover the analog of the electric 
field is the Newtonian or gravitoelectric field g, while the 
analog of the magnetic field B = V x A is the gravit¬ 
omagnetic (or “frame dragging”) field Q, both defined 

by 

g = -V0, Q = V x h. (2.10) 


For many systems of interest the GEM equations are 
easily solved in the same way as those of EM, for example 
by the use of a Biot-Savart type law. For a stationary 
spinning sphere like the earth or (a ppr oximately) a neu¬ 
tron star the exterior solutions are j20| i0 ' 



r 2 GI _ „ 

h = ——lo x r, 



n = 


2 GI 


[3(r • w)r — w], 


( 2 . 11 ) 


where u is the spin and / is the moment of inertia of the 
sphere; for a uniform density sphere / = (3/5 )Mr 2 s 0- 
The fields in Eq. (2.11) are of course time independent. 

Another body of theoretical interest is a hollow spher¬ 
ical shell, since the gravitoelectric field vanishes in its 
interior, while the gravitomagnetic field near the center 
is about 0 


n = 


2GM _ 

- ui. 

Ths 


( 2 . 12 ) 


A gyroscope at the center of such a shell would precess 
at H/2 0. 


3. SCHROEDINGER EQUATION FORM (SEF) 
AND ELECTROMAGNETIC INTERACTION 


Equations (2.6a) are the field equations, while Eq.s (2.6b) 
follow from the Lorentz conditions. For low velocity 
sources we expect 4> ~ <p(v/r) and <fi ~ <p(v 2 /r 2 ) where r 
is the characteristic distance to the source, so we will ig¬ 
nore (j> in Eq.s (2.6) and everywhere else henceforth, but 
we will retain <f>. Similarly we will ignore both h° and h J . 
Then Eq.s (2.6) become 

V 2 ^ = 4t rGp, V' 2 h J = -16t rGpv J , 40 - V • h = 0. 

(2.7) 

We refer to Eq.s (2.7), with slowly varying fields, as 
the GEM (for gravitoelectricmagnetic) approximation or 
limit. 


First we briefly discuss the free scalar Klein-Gordon 
equation and recast it into a Schroedinger equation form 
(SEF), which is useful in obtaining the non-relativistic 
limit 00. The SEF is exact and equivalent to the 
usual Klein-Gordon equation. Surprisingly, it does not 
appear to be well-known and we have found no reference 
to it in the literature 00. 

The free Klein-Gordon Lagrangian and equation are 

L = (p — V 2 <£ + m 2 ip = 0. (3.1) 

The time variation of the scalar field ip due to the rest 
energy m can be separated out by the transformation 

ip = e ~ imt ip. (3.2) 
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Substitution of this into Eq. (3.1) leads to an equivalent 
form for the Klein-Gordon equation, 


. • V 2 ip ip 

zip = -1-. 

2m 2m 


(3.3) 


Equation (3.3) is exact and has the form of a non- 
relativistic Schroedinger equation, but contains an ex¬ 
tra term, a second derivative, hence our appellation of 
SEF. An exact solution to Eq. (3.3) is the plane wave 
exp(—iE n t + ip ■ x), where E n = E — m is the non- 
relativistic energy. 

In a low velocity system the time variation of ip is asso¬ 
ciated with non-relativistic kinetic energy,which is much 
smaller than that of e _lmt associated with rest energy. 
Thus the last term in Eq. (3.3) will be small for posi¬ 
tive energy solutions, and can be viewed as a relativistic 
correction; for negative energy solutions this is obviously 
not so. To lowest order approximation the relativistic 
correction corresponds to the standard p 4 term in the 
expansion of the relativistic energy. 

Note that if the second time derivative in Eq. (3.3) is 
treated as a small perturbation the Cauchy initial value 
structure is that of a first order Schroedinger equation. 
Of course this presents no contradiction with the Cauchy 
structure of the second order Klein-Gordon equation - 
so long as we limit ourselves to positive energy solutions 
0- 

To include electromagnetic interactions we use the 
standard gauge in variant minimal coupling recipe, id M —>■ 
13( 241. The Lagrangian then becomes 


idy, - eAy 


L = - eAy<p*){iip\ v - eA v ip)rf v - mVV, 

(3.4) 

which gives the Klein-Gordon equation in manifestly 
gauge invariant form, 


(idy - eAy)(id v - eA I/ )Tf u <p - m 2 ip = 0. 


(3.5) 


Equivalently we could apply minimal coupling to the 
free SEF Eq. (3.3) to obtain Eq. (3.7). Eq. (3.7) is 
an SEF for the spin zero particle interacting with the 
electromagnetic field in the standard way, with a second 
derivative relativistic correction. It is exact and useful in 
obtaining the non-relativistic limit, but apparently does 
not appear in the literature [22j. Note that in Eq. (3.7) 
there occurs a term ieS7 Aip in the kinetic energy and a 
term iV ip in the relativistic correction. In the Lorentz 
gauge that we are using they cancel and are not really 
present in the SEF of Eq. (3.7). 

Eq. (3.7) contains a somewhat subtle feature related 
to conservation of probability and reality of the energy, 
which we next consider. In the Schrodinger equation of 
nonrelativistic quantum mechanics, 

p 2 

iip = Hip = ——ip + Vip, (3-8) 

2m 

the potential V is usually taken to be real, thus mak¬ 
ing the Hamiltonian Hermitian; this guarantees that its 
energy eigenvalues are real and that probability is con¬ 
served. An equivalent way to analyze probability conser¬ 
vation is in terms of a probability density and associated 
3-vector current. These are defined as 

P = 3 = ^ (ip*Vip - . (3.9) 

From Eq. (3.8), we readily obtain the conservation equa¬ 
tion 


p + V-j=iip*(V-V*)ip = 0. (3.10) 

Thus probability is conserved if the potential is real. 

For the relativistic Klein-Gordon equation interacting 
with the electromagnetic field the analysis of probabil¬ 
ity conservation proceeds in a similar way. The 4-vector 
current is defined as 


Equation (3.5) may also be written in a form displaying 
the electromagnetic coupling separately on the right side, 

V>\n\vrf v + m 2 ip = - 2ie<£| M A,r/ i!/ - ieA il \ v rf lv ip 

+ e 2 Ay,A v r)^ v ip. (3.6) 


Notice that we use only lower indices on the potential 
A v and only upper indices on the Lorentz metric r) tlv , 
which will prove to be convenient. Also for convenience 
we use the Lorentz gauge, without loss of generality, since 
the Lagrangian is gauge invariant. Thus henceforth we 
take = V + eV • A = 0, so the second term 

on the right of Eq. (3.6) is zero. The Lorentz gauge 
is particularly appropriate for problems involving a time 
independent potential or an external radiation held. 

To recast Eq. (3.6) into SEF we make the substitution 
Eq. (3.2), as before, to obtain 


i^=^L+Vif>- ^ ,y + eV- A = 0 (3.7) 


2m 
V = eA°, 


2m 

p = *V, n = p — eA. 








\ eA^ 

1 - 

/ m 


p <p 


i eA^ 

= —p d V'ip - p p. 

2 ?7i m 


(3.11) 


Conservation follows from the Klein-Gordon Eq. (3.5) or 
(3.6), and is expressed as ji'i = 0. The same j M , up to a 
constant factor, follows from the Lagrangian Eq. (3.4) as 
the electromagnetic current. We normalize the current so 
that the 3-vector part is consistent with the Schrodinger 
current. In terms of ip defined in Eq. (3.2) the density 
and 3-vector current are 


3° = ( 1 - —) 

\ my 2m 

n ^p A 

j k =—iP* d ^ip - - — ip*ip. (3.12) 

2m m 


Thus the conserved probability density j° is not the same 
as the Schrodinger density ip*ip in Eq. (3.9), although it 
is simply related to it, as we now discuss. 


7feol4e«taiZc.a/ P/ty.lZc . 1 









4 


It is instructive to consider the case of an energy eigen- 
funtion ip(x,t) —> e~ lEt ip(x). Then the probability den¬ 
sity in Eq. (3.12) becomes, approximately, 


f = 1 + 


E —V 


%p*{x)ip(x) = ( 1 + 


2 TO 2 


ip*(x)lp(x). 

(3.13) 


We can view Eq. (3.13) as a relativistic renormalization 
of the wave function since the factor (1 + p 2 /2m 2 ) is a 
Lorentz volume contraction factor. This is easy to see if 
we write it in terms of the velocity, defined as v = p/m , 


correct to order 0(mv 4 ), 


Eips =(1 +p 2 /4m 2 ) 

„4 


(lm + V ) t 1 _ 


P ^ s+ 2 to [^ ,p2 / 2m ] ^ s 


8m 3 

J2 


= (jj^ +V^jip s + [p 2 / 4m 2 , V] ip s - 


8m 3 


ips 


+ ^ [F,p 2 /2m] ip s 


p 




8m 3 


ips- 


(3.18) 


-2 -*2 
1 + 2^~ 1 + y = 


x/T^ 


(3.14) 


We will illustrate the importance of the wave func¬ 
tion renormalization for the important example of a 
slowly moving particle in a time independent potential 
eA° = V and zero 3-vector potential - that is V = 0 
and A = 0. This calculation is the spinless analog of the 
well-known hydrogen fine structure calculation based on 
a non-relativistic reduction of the Dirac equation[25|. To 
seek energy eigenvalues we again set ip{x, t ) —> e~ lEt ip(x ) 
and Eq. (3.7) becomes 


Thus the non-Hermitian term has been cancelled by the 
wave function renormalization, so the resulting equation 
will have real energy eigenvalues. It is worth noting that 
the procedure leading to Eq. (3.18) is analogous to that 
leading to the hydrogen fine structure energy, but the 
result is simpler since there are no spin-orbit and Darwin 
terms [25j. 

In summary, Eq. (3.7) is gauge invariant, and contains 
only Hermitian operators when the gauge is carefully cho¬ 
sen and wave function renormalization is applied. Similar 
considerations will be used in sections 4 to 6. 


4. GRAVITATIONAL INTERACTION 


(E-V)iP = ^-iP--^-(E-V) 2 ip. (3.15) 

2m Am 

The operators on the left side and the first term on the 
right side are of order 0(mv 2 ) and the last term on the 
right is 0(mv 4 ); we will work to only this order. We 
iterate Eq. (3.15) and obtain 


E^=^ + V^-±(E-V)A- r 


~*P + Vip — [(E — V),p 2 /2m] ip 


1 


4?7l 2 


(E - V)ip 




ip + ^ [V,p 2 /2m] ip. (3.16) 


However the last term in Eq. (3.16) is not Hermitian, 
as is easily verified, so it is not a consistent eigenvalue 
equation for the real energy. This is because the wave 
function ip is not correctly normalized. We renormalize 
it according to Eq. (3.13) by defining a wave function 
that is normalized like the Schrodinger wave function, 

ips — (1 +p 2 /4m 2 )i/’, ip = (1 - p 2 /Am 2 )ips- (3.17) 


Substituting Eq. (3.17) into Eq. (3.16) we obtain for ipg, 


For a scalar field the interaction with gravity is ob¬ 
tained by making the Klein-Gordon equation generally 
covariant. To do this we need only multiply the flat space 
Lagrangian by the square root of the absolute value of 
the metric determinant (denoted ^/g) to form a scalar 
density, replace the Lorentz metric by the Rieman- 
nian metric g ^, and substitute covariant derivatives for 
ordinary derivatives [26j . Then, in covariant form, the 
Lagrangian in Eq. (3.4) becomes a density, 

L = - eA^p*){ip\ v - eApip) g^" - mVV], 

(4.1) 


and the Klein-Gordon Eq. (3.6) becomes, 

<PMv& lv + ™V = - Ziep^Avg 11 '' - ieA^g^p 

+ e 2 A„A v g^ip. (4.2) 

The double slash here denotes a covariant derivative and 
the single slash an ordinary derivative; for a scalar they 
are the same. The covariant dAlembertian and diver¬ 
gence terms are conveniently expressed as fisj ]. 

(y/g<P\»gr)\ v , 

A^ = A^ v =-L(y/gA v )\ v . (4.3) 

We have here assumed only the simplest couplin g t o 
gravity, consistent with the equivalence principle [211 ]. 
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Additional terms involving the Riemann tensor are also 
possible since they vanish in flat space and are thus con¬ 
sistent with the flat space Klein-Gordon equation 0 • 

Covariant electromagnetism is gauge invariant so we 
again choose for convenience the Lorentz gauge, Ay = 
0; the second term on the right of Eq. (4.2) is then 
zero. Recall that we also using the Lorentz gauge for the 
gravitational field, Eq. (2.4). 

We next substitute Eq. (4.3) into Eq. (4.2), expand 
the metric to first order, g^ 1 ' = r]^ — , and apply the 

Lorentz gauge condition Eq. (2.4); this yields 

v"if v +rri 2 <p 

=1 P\ij.\v^ u + (—+ e 2 A^Avip)^^' — /d" ). 

(4.4) 

This surprisingly simple result Eq. (4.4) only holds in the 
Lorentz gauge for both the electromagnetic and gravita¬ 
tional fields. Eq. (4.4) holds to all orders in the elec¬ 
tromagnetic field, and to first order in the gravitational 
field. The right side includes the electromagnetic terms of 

(3.6) , an additional term for coupling to the gravitational 
field, and two mixed terms for coupling to the product of 
the electromagnetic and gravitational fields. The mixed 
terms are novel and may be of possible interest in con¬ 
nection with very strong electromagnetic fields, such as 
those found in laboratory magnets or intense laser light, 
as we will discuss in section 6. 

To recast Eq. (4.4) in SEF we separate out the rest en¬ 
ergy time dependence using Eq. (3.2) as before. Straight¬ 
forward manipulation leads to a generalization of Eq. 

(3.7) to include gravitational interaction terms, which we 
write as 

+ + + (4.5) 

The term L/g0 contains interactions with only the gravi¬ 
tational field, and in the GEM approximation of section 
2 it is 


Eq.s (4.5) to (4.7) are a restatement of Eq. (4.4) in terms 
of the wave function 0 and in the GEM approximation for 
the gravitational field. The only approximations made 
so far are the GEM of section 2 for the external gravi¬ 
tational field, which assumes a slow moving source. The 
electromagnetic field is treated exactly and there is no 
low velocity approximation for the particle. 


5. NON-RELATIVISTIC LIMIT 


In this section we consider the low velocity limit of 
Eq.s (4.5) to (4.7) to obtain our main, remarkably simple, 
result. The four terms in the first parenthesis of Eq. 
(4.6) are 0 times, respectively, 0(m), 0(mv 2 ), 0(mv 2 ), 
and 0(mv 4 ); we will work to 0(mv 2 ) and drop the last 
term. In the second parenthesis the terms are h times, 
respectively, 0(mv) and 0(mv 3 ); again we drop the last 
term. We then do the same order counting in Eq. (4.7) 
and combine the remaining terms with those of Eq. (4.6) 
to obtain an SEF, 


*0 = 


n 2 0 

2 m 


0771 


+ E0- 


(it% - E) 2 0 


2 m 


2n 2 \ , 

! + —r 0- 


m* 


h • ffy. 


(5.1) 


Thus the Newtonian field 0 couples to the particle mass 
with a small correction factor, and the gravitomagnetic 
field h couples to the kinematic momentun n. We will 
discuss the gravitomagnetic coupling at length in the 
next section. 

The factor preceding the Newtonian field in Eq. (5.1) 
is interesting. It may be written in terms of the particle 
velocity as (1 + 2v 2 ), with the velocity defined as v = 
n/m. We can show that (l+2v 2 )0 is simply the potential 
that the moving particle sees in its rest frame. To see this 
we apply a Lorentz transformation in the x direction to 
the metric perturbation, 


H g il> = 0 f m0 - — + 2 (*0) - ^ 

\ rn m 


+(h-p) 0 + 


(iip) 


(4.6) 


The mixed term U ge 0 contains interactions with the 
product of the gravitational and electromagnetic fields, 
and in the same approximation it is 


7W=(- 2 W -^)-20j4lW) 

V m m 

m m / 

( j V 2 ip e 2 A 2 ip Ve(A-h)ip\ 

+ 0-b0-1-• 

\ m m ml 


(4.7) 


ti^ = a» a a v ph a0 , 


7 = 


y /1 — V 2 


( 7 -V7^\ 

\ —v 7 7 ) 


(5.2) 


Ignoring the small off-diagonal terms of the metric we see 
from Eq. (2.5) that h 00 and thus the Newtonian potential 
transform like 


0' = (a° o a° o )0 + (a 0 ! 0^)0 = q 2 (l + v 2 )0 “ (1 + 2v 2 )0. 

(5.3) 

Finally, we emphasize that Eq. (5.1) contains only 
Hermitian terms in the gravitational interactions; recall 
also that in section 3 we showed that the rest of the 
equation contains only manifestly Hermitian terms if the 
Lorentz gauge is used and the wave function is properly 
renormalized; thus no problem arises with conservation 
of probability and reality of the energy eigenvalues. 
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6. GRAVITOMAGNETIC PHYSICAL EFFECTS Eq.s (5.1) and (6.1), 


We will briefly discuss some physical effects implied 
by Eq. (5.1), concentrating on gravitomagnetism. Our 
purpose is only to provide some simple parametrizations 
and sample numerical estimates. A specific and detailed 
discussion of any laboratory experiments or astrophysical 
observations is beyond our present scope, so we leave it 
as a challenge to experimentalists to design and analyze 
specific experiments. 

We will consider here only neutral systems, perhaps 
containing charged particles, since electromagnetic ef¬ 
fects on a system with net charge generally swamp grav¬ 
itational forces by many orders of magnitude [I]. 

The gravitoelectric or Newtonian term (prmp in Eq. 
(5.1) is generally much the dominant gravitational in¬ 
teraction. Its effects have been seen in diverse experi¬ 
ments, notably on neutrons in the earth’s gravitational 
field, which make it clear that nucleons do indeed fall and 
behave as expected [8j][§]! That term is also responsible 
for most of the phase shift seen in recent atomic beam in¬ 
terferometry gravity experiments 0- Such experiments 
may be able to test the equivalence principle to 10~ 15 or 
better; however the phase shifts, including small subtle 
corrections, may be calculated more easily and accurately 
with the semi-classical approach of Dimopoulous et. al. 
that treats the atoms as neutral particles, rather than 
the approach we adopt here 0 . 

Our present interest is more on the gravitomagnetic 
term h ■ II ip of Eq. (5.1). The first part, h ■ pip, is the 
analog of the EM coupling to the vector potential, as we 
see by writing out the kinetic energy term [27}, 

ft 2 ip fPip eA -pip + ieS7 ■ At)) + e 2 A 2 ip 

2m 2m m 2m 2m 

From the second term of Eq. (6.1) the correspondence 
between magnetic and gravitomagnetic couplings is thus 
evident, 


r e y 

h <^>- A. 

m 


( 6 . 2 ) 


It is illuminating to take the quantum system to be a 
fictitious atom containing a scalar electron of mass m e , 
which we then place in B and/or f l fields that are ap¬ 
proximately constant over the size of the atom. In this 
case we may conveniently choose the potentials as 


g -* g —, g 

- A ■ pip = —— —B x r ■ pip = — ——B ■ r x pip 

m e 2m e 2 m e 


2m, 


■B ■ Lip , 


h ■ jhp = — Q x r • pip = — O • r x pip = -Al ■ Lip. 


(6.4a) 

(6.4b) 


Here the orbital angular momentum L is taken with re¬ 
spect to the center of mass. Thus Eq. (6.4a) means 
that orbital angular momentum produces a magnetic mo¬ 
ment ft = (e/2 m e )L, since Eq. (6.4a) is the energy of 
a magnetic moment in a magnetic field. In the same 
manner we interpret Eq. (6.4b) to say that orbital an¬ 
gular momentum produces a gravitomagnetic moment, 
jlgrav = (— 1/2)L. It is clear from the correspondence 
that, since a magnetic moment precesses at the Lar- 
mor frequency eB/2m e in a B field, the gravitomagnetic 
moment will precess in a gravitomagnetic field with fre¬ 
quency f2/2, and in the opposite direction. 

The classical analog of the gravitomagnetic coupling in 
(6.4b) leads to the Lense-Thirring precession of a gyro¬ 
scope in the field of the earth as observed in the Grav¬ 
ity Probe B experiment @0- Thus a quantum system 
such as our fictitious atom should undergo the same pre¬ 
cession as a classical system in a gravitomagnetic field. 
This should be expected for a fundamentally geometric 
effect, according to general relativity. 

We note that a similar analogy holds for classical or¬ 
bital motion: a charged particle in a constant B field 
(and no E field) orbits at the cyclotron frequency eH/m e , 
whereas a massive particle in a constant f2 field (and no 
g field) orbits at frequency f l [2l[. It is of course easy to 
produce a region of B = constant E = 0 in the lab; it is 
similarly easy to produce a region of = constant g = 0 
by using the interior of a massive spherical shell, as we 
noted in section 2 [2lJ. Thus a particle near the center 
of such a hollow shell will orbit at H = (:2GM/r/, s )w. 

It is of some interest to estimate the gravitomagnetic 
precession frequency of the fictitious atom near the sur¬ 
face of two bodies of interest, the earth and a neutron 
star, both approximated by a spinning sphere with a 
gravitomagnetic field given by Eq. (2.11). In order of 
magnitude 


n 


2G(Mrl) 


' S' b , 


(6.5) 


A = —B x r, h = -0 x r, (6-3) 

where f = 0 is the center of mass of the atom. The 
magnetic and gravitomagnetic interactions are then, from 


where r s = 2GM is the body’s Schwarzschild radius and 
rt, is its actual radius. For the earth, r « r;, = 6 x 10 3 km, 
r s = 1.8cm, oj = 7.3 x 10~ 5 rad/s. For a typical rapidly 
spinning neutron star we take r ps rb ~ 10km, r s « 
lkm, ui = 10 2 rad/s. Then the precession frequencies and 
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associated energies are of order 

12 « 10 _13 rad/s, ifchar ~ M2 « 10~ 28 eV (earth ) 

(6.6a) 

fl 10rad/s, ifchar ~ M2 « 10 _14 eV (neutron star) 

(6.6b) 


The energies in Eqs. (6.6) are of course extremely small 
compared to typical atomic energies, even hyperhne ener¬ 
gies; in any situation where GEM effects directly compete 
with EM effects they will naturally be swamped by many 
orders of magnitude. 

The other part of the gravitomagnetic term in Eq. 
(5.1), -h-eAip, might conceivably be of interest for some 
terrestrial or astrophysical environments if the vector po¬ 
tential A is sufficiently large. Two example systems of 
interest, obtainable in a terrestrial laboratory, are a high 
field magnet and an intense laser beam [28[ [29| . We first 
consider a quantum system such as our fictitious atom in 
a large B field. The perturbation energy, with the use of 
Eq. (2.10), is of order 


A E « (h ■ eA) « f2e.Br 2 



(6.7) 


Here r q is the characteristic size of the quantum system, 
A c = H/m e c « 2.43 x 10~ 12 m is the Compton wavelength 
of the electron, and we have included a factor of c to 
make the dimensions clear. For a numerical estimate we 
take r q w 10 _9 m, 12 « 10~ 13 rad/s from Eq. (6.6a), B « 
10 2 T, and express the Bohr magneton as e/2 m e = (5.8 x 
10 _5 eV/T). Then the energy and associated frequency 
are of order, 

A E « 10 _32 eV, Aw « 10 _17 rad/s (atom in terrestrial B) 

( 6 . 8 ) 


This is even smaller than the values in Eq. (6.6a) so 
the presence of the magnetic field would not appear to 
produce a larger effect of interest. 

However if we suppose a much larger quantum system, 
perhaps Cooper pairs in a Josephson junction of nearly 
macroscopic size r q « 10 _6 m, then the above values in¬ 
crease by six orders of magnitude and possibly become 
more interesting [30], 

A E « 10 _26 eV,Aw « Hr n rad/s 

(large quantum system in terrestrial B) 

(6.9) 

We leave it to experimentalists to consider such larger 
systems. 

Finally we consider a fictitious atom in an intense laser 
beam of frequency w; and electric held E. For a radiation 
Held A ss E/uji and we may estimate the perturbation as 


before 

A Ek, (h- eA) 

= ( eEl ' 
\m e iOi 




n. 


( 6 . 10 ) 


The quantity eEi/uJi is the approximate energy trans¬ 
fer to the electron in one cycle of the laser beam, so 
the dimensionless parameter ag = eE/m e u>i is the frac¬ 
tion of the electron rest energy absorbed from the laser 
beam in one cycle; it is a measure of how relativistic 
the quantum system can become, and is thus also a con¬ 
venient laser intensity parameter [2f|. Note that we 
may also express ao in terms of the Schwinger critical 
Held, at which spontaneous pair production is relevant, 
as ao = ( E/E S c){m e /ui ) [31j. 

If we take as rough examples ao = 1 and r q « 10~ 9 m 
we obtain from Eq. (6.10) 

Aw « 10~ n rad/s , A E « 10 _26 eV 

(atom in terrestrial laser) (6.11) 


Of course our approximations do not hold for ao much 
larger than 1, and are probably not even very accurate 
for ao of order 1. The values in Eq. (6.11) are larger 
than those in Eq. (6.6a) so the presence of the laser field 
may increase the gravitomagnetic effect and possibly be 
of interest to experimentalists. 

For a larger quantum system, for example a Josephson 
junction as mentioned above, the values might increase 
by perhaps three orders of magnitude, according to Eq. 
(6.10), to about 

Aw « 10 _8 rad/s , A E « 10“ 23 eV 

(larger system in terrestrial laser) 

( 6 . 12 ) 

We note incidentally that measurements of the Joseph¬ 
son constant 2 e/h (relating vo ltag e to frequency) reach 
a precision of about 2 x 10~ n [331 ]. 

There is a major difficulty associated with detecting 
gravitomagnetic effects that we have not yet mentioned; 
any rotation of the laboratory apparatus will in general 
compete with and swamp gravitomagnetic effects via the 
Sagnac effect [34j]. Thus any experiment would need to 
provide a way to suppress or separate such rotational 
effects at about the level of 10 -13 rad/s according to Eq. 
(6.6a). 


7. SUMMARY AND CONCLUSIONS 

Our study of the interaction of a spinless quantum 
particle with an electromagnetic field and a weak gravi¬ 
tational field has lead to a rather simple non-relativistic 
limit expression that is useful for considering possibly ob¬ 
servable gravitomagnetic effects. The laboratory obser¬ 
vation of such effects is clearly a very difficult prospect, 
as our illustrative examples indicate, but would be of 
fundamental interest. 
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